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Vortices with fractional flux 
in two-gap superconductors and in extended Faddeev model 
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We discuss linear topological defects allowed in two-gap superconductors and equivalent extended 
Faddeev model. We show that in these systems there exist vortices which carry an arbitrary fraction 
of magnetic flux quantum. Besides that we discuss topological defects which do not carry magnetic 
flux and describe features of ordinary one-magnetic-flux-quantum vortices in the two-gap system. 
The results could be relevant for the newly discovered two-band superconductor MgB2- 



A fundamental property of the Abelian Higgs model 
is the quantization of magnetic flux 0|. In an ordinary 
superconductor the Abrikosov vortices can carry only in- 
teger number of magnetic flux quanta. The intriguing 
possibilities of topological defects carrying a fraction of 
flux quantum have long attracted interest and several 
nontrivial realizations were identified. For example, a 
half flux-quantum vortex in a spin-1 condensate is a con- 
figuration where a Cooper pair has its spin reversed when 
moving around the vortex core (this is an analogue of an 
Alice string in high energy physics where a particle mov- 
ing around the string flips its charge or enters a "shadow 
world" ) also a half fiux-quantum vortex can be formed 
on a junction of 3 grain boundaries in a crystal etc In 
this paper we discuss vortices in two-gap superconductors 
|3,|j (known in particle physics as a Higgs doublet model 
|5|) and in the extended Faddeev model. We show that 
these vortices can carry an arbitrary fraction of magnetic 
fiux quantum. 

Experimentally, two-gap superconductivity has been 
observed in the transition metals Nh, Ta,V and in Nb- 
doped SrTiOz |^ . More recent experiments indicate the 
two-gap nature of superconductivity in MgB2 and 
2H — NbSe2 ^ ■ Two-gap models appear also in the the- 
oretical studies of liquid metallic hydrogen, which should 
allow superconductivity of both electronic and protonic 
Cooper pairs |^ . In liquid metallic deuterium a deuteron 
superfluidity may be present along with superconductiv- 
ity of electronic Cooper pairs Q. Other realizations of 
the two-gap system are superconductors with two types 
of pairing (e.g. a mixture of s- and p-wave condensates). 

A two-gap superconductor is described by a two- 
flavour Ginzburg-Landau free energy functional: 

F= |(V-f ieA)*i|2 ' ^ 
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1,21 



2 

|2^ 



ieA) ^'2 



(1) 



^I^'qI* and 77 is a characteristic of the interband Joseph- 
son coupling strength Q . 

Many exotic properties of (|l|) are obscured in the 
Ginzburg-Landau presentation of the free energy func- 
tional. In fioll it was shown that there exists an exact 
equivalence mapping between the model (^) and an ex- 
tended version of Faddeev's 0(3) nonlinear cr-model [ pr] , 
which describes the two-gap superconductors in terms 
of gauge invariant variables which explicitly show the 
degrees of freedom present in the system. This model 
consists of a three-component unit vector n in interac- 
tion with a massive vector field C and a density-related 
variable p [n2|: 



F = ^(Vn)2 + [Vpf + + V{p, n3) -f p^Kn^ 
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wheied, = = A + Bns- 
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-Cn|. The models @ and 



(g) are connected in the following way ||1C 



&2'7i2]; B ~ p^[8c2m| — 8ciml — 52m2 



coefficients 

A, B, C are given by: A — p^[Acimf + 4c2m2 — 5imi — 
p^[8c2m| — Scimf — 52m2 -I- 6imi]; C — 
The position of the unit vector fi on 
the sphere S'^ can be characterized by two angles as fol- 
lows: n = (sin 6* cos 7, sin sin 7, cos 0), where 7 = (^i — 
02);|^'i,2| [\/2^psin(f) ,V27^pcos(f)]; C = 
^ {M-tV*! - vPiVM-t} + ^{vE'*VvI/2 - *2VVI/*} - 
2| _|_ \_^2\_\ consider the system in the 

p \ mi m2 y ^ 

presence of an interband Josephson coupling 77[\E'*\I'2 -I- 
^2*1] p^Kui where K = 2?7^toiTO2- The poten- 
tial term V in (||) determines the energetically preferred 
ground state value for rt3, which corresponds to uni- 
form density of both condensates. We denote it as 

r 1 r 1 ~i 

^3=^-^ iVi_^jV2 ^ stands for the 
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average < |^'i,2p >■ Thus, in the presence of the intrin- 
sic Josephson effect, the ground state value of n corre- 
sponds to a point where ni = —1 on a circle specified 
by the condition 77.3 = 713 = cos 6'. We note that (||) 
shows the Meissner effect, which is manifest in genera- 
tion of the mass for C The corresponding length 
scale A is the London magnetic field penetration length: 



+ 



1*2 



. One can observe from (0) that 
along with the Meissner effect in a two-band supercon- 
ductor there exists a neutral boson. That is, if one con- 
siders a uniform density of condensates, the vector field 
C decouples from the field n, because n • 9^11 x djii cx 
sni9[di'^dj9 — di^djO] — 0, when 713 = cos^ = const. 
Thus, when 77 = 0, the system possesses a massless neu- 
tral 0(2) excitation associated with the phase variable 
^ — {(pi — 4>2) ■ This phenomenon has no counterpart in 
one-gap superconductors where the Goldstone boson as- 
sociated with 0(2) symmetry is "eaten" by the Meissner 
effect. In the presence of the Josephson effect (rj ^ 0) 
the variables still decouple but the Josephson term p^Kni 
breaks the neutral 0{2) symmetry by giving a nonzero 
mass to ni. 

Let us discuss vortices in the system (^. In the Lon- 
don hmit (l^'il = const), the eqs. (|l|), (0) become: 
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sir? e{V-if + p2 ^sin^l^- j V01 -I- cos^ j V(/)2 - eA 
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H — Y + sin 9 cos 7 



(3) 



From this expression one can observe that the vor- 
tices characterized by A((/)i -I- (f)2) = §^ d/[V((/)i -I- ^2)] = 
Airm; A((/)i — (j)2) — 0, (where we integrate over a closed 
curve a around the vortex core) is the analogue of tti-Aux 
quanta Abrikosov vortices in an ordinary superconductor 



characterized by 



1*1^ _ 



\ mi 



Let us observe 



that if both phases ^1,2 change by 2tt around the core 
then a vortex carries one quantum of magnetic flux. 

The vortices characterized by A(0i — (t>2) = 4?™ in 
the case where ^'^^ ^ l*iL have nontrivial structure. 

mi ' m2 

Let us first consider the case when A(0i -I- 4)2) = 
and A (01 — 02) = 47r. First of all such a vortex fea- 
tures a neutral superfiow characterized by a 47r gain 
in the variable 7. In the case of 77 7^ this vor- 
tex is described by the sine-Gordon functional: F 



(l/4)p2 sin^ 9{V^Y+p^K sin 9 cos 7. [In the hmit 

1*2 



Ifil" 



-, this vortex does not involve a nontrivial config- 
uration of the field C, as can be seen from (^)]. In 
the case = 77 = this vortex is equivalent to a 
vortex in a neutral superfluid with superfluid stiffness 
(l/2)/9^ sin^ 0. As follows from (Qj, in the case when 



associated with 7 is necessarily accompanied by a non- 
trivial configuration of the charged field C. That is, for 
such a topological defect the second term in (^ becomes 

v(0i-fe) „,a1" This term 



{sin^(|)-cos^(f)} 



P 

is nonvanishing for such a vortex configuration, which 
means that this vortex besides neutral vorticity also 
carries magnetic field. Lets us calculate the mag- 
netic flux carried by such a vortex. The supercurrent 
around the core of this vortex is: J = 2ep^[{ sin^ (|) — 
cos'^ (I) } ^^'^^~'l"^'> — eA] . Lets us now integrate this ex- 
pression over a closed path a situated at a distance larger 
than A from the vortex core. Indeed at a distance much 
larger than the penetration length the supercurrent J, or 
equivalently the massive fleld C, vanishes. Thus we arrive 
at the following equation: 



<i> — cos ( 



^V(0i 



62)fil = cos 0$o = n-3$o, (4) 



where ^ = §^ Adl is the magnetic flux carried by our 
vortex and $0 = 27r/e is the standard magnetic flux 
quantum. From (Qj it follows that in our system such a 
vortex can carry arbitrary fraction of magnetic flux quan- 
tum since it depends on the free parameter cos 9 which is 
a measure of the relative densities of the two condensates 
in the system. In the general vortex charac- 

terized by the following phase changes around the core 
A01 ~ 27rfci and A02 — 27rfc2, carries the following flux: 



(kiM) 



Qki+C0S^Qk2 



(5) 



, a topological defect in the neutral field 



We should emphasis that a straightforward inspection 
of (H) shows that this model allows neutral vortices asso- 
ciated with the neutral 0(2) boson without a nontrivial 
configuration of the field C for any values of 77.3. However 
such solutions (e.g. vortices in the case 713 ^ charac- 
terized by A7 = 27r; C = 0) are unphysical because these 
vortices do not satisfy the condition that 0^ change by 
2Trki around the vortex core (as follows from (||)). Thus, 
while the original model (|^) and the extended Faddeev 
model (^) have the same number of degrees of freedom 
so that the fields n and C are dynamically independed by 
construction, the mapping incurs a constraint on topo- 
logical defects in n and C since 0i_2 appear in both of 
them. So, in the hydrodynamic limit in (^ there is no 
direct coupling between the fields n and C, however in 
a non-simply-connected space a topological defect in the 
field n necessary induce a nontrivial configuration of the 
field C. The consequence of this is the fractionalization 
of magnetic flux in the model (^). 

So a two-band superconductor (|l|) and extended Fad- 
deev model (||) allow the following linear composite topo- 
logical defects: 

(i) A(0i — (j>2) = 47r77; A(0i -I- ^2) = 0: these are the 
vortices which feature neutral superflow. When the 
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ground state of the variable n corresponds to the equator 
on (that is the case when i^ii- = I^lL or equivalently 

^ mi 7712 

cos = 0) these vortices do not carry magnetic flux. In 
the case when cos 0, the only physical solutions with 
neutral vorticity also carry a fraction of magnetic flux 
quantum. In the case of nonzero Josephson coupling 
these vortices are described by the sine-Gordon equation 



(ii) A7 = 0; A(0i 



= Ank these are the vortices 



which feature circular supercurrent and no circular neu- 
tral superflow. Such vortices carry k flux quanta of mag- 
netic field. One of the physical consequences of our 
analysis is that we can observe that in the two-gap super- 
conductor in spite of the existence of two types 
of Cooper pairs, a formation of two sublattices of vortices 
corresponding to each of the condensates in an external 
field is energetically forbidden. This is because the en- 
ergy per unit length of noncomposite vortices is divergent 
in an infinite sample both in cases of zero and nonzero 
Josephson coupling (in case of finite 77 a vortex creates a 
domain wall which makes its energy per unit length diver- 
gent in infinite sample mjl^/ j. We can also observe that 
a superconductor made up of one condensate which is of 
type-II and another condensate which is of type-I will 
in general, in the presense of an external magnetic field, 
form one flux quantum vortices involving both conden- 
sates and will preserve two-gap superconductivity, even 
if the external field exceeds the thermodynamic critical 
magnetic field for the type-I condensate. 

(iii) A7 = 27rn; A{(f>i -\- <j>2) = 47r^: These are the vor- 
tices which may be viewed as a co-centered Z-flux quanta 
Abrikosov vortex combined with a vortex with neutral 
vorticity with winding number n carrying a fractional 
magnetic flux. These vortices and the vortices of type 
(ii) characterized by A{(j>i -I- (^'2) > Stt could be unstable 
against decay into more simple vortices [ p^ . 

Besides composite vortices the system allows vortices 
where the phase of one of the condensates changes by 
27T around the core whereas the phase of the second con- 
densate remains constant. These vortices are also prin- 
cipally different from ordinary Abrikosov vortices. Let 
us discuss these vortices in the Ginzburg-Landau formal- 
ism (Q). This will provide additional perspective on the 
discussion and will allow us to illustrate the connection 
between the models (|l|) and (||) . We shall present a solu- 
tion for these vortices in the limit when the inverse mass 
for Til which is given by Josephson term is larger than 
the sample size. Such an approximation allows one to 
put 77 to zero. We should note that there certainly ex- 
ist systems where 77 is exactly zero. That is, a model 
with 77 = should describe vortices in a bi-layer system: 
a superconductor-insulator-superconductor compound is 
an example of a system of two condensates coupled only 
by a gauge field. In such a system, a vortex in one layer 
carries flux through the second layer, where it also in- 



duces a current, which should lead to the flux fraction- 
alization. Besides that an additional U{1) symmetry ap- 
pears if one considers equal-spin pairing. In that case, if 
the spin-orbit coupling is neglected, each spin population 
has its own phase, without Josephson coupling | p5| . 

Let us recall the GL equation for the supercurrent 
in standard notations: J = 7^[4'JV5'i - ^'iV^^] -f 



nation when only the phase 
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. In the sit- 
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1 changes by 27r around 
the core while the phase of the second condensate re- 
mains constant, we have the following expression for 

the vector potential: A - 1 ^2 i*iP i*2jf 
-1 



1 i*ir 

e mi 



i*ir 



-t- 



i*2i- 



7771 ^2 

lows that the vortex characterized by A0i 
carries the following fractional magnetic flux 



V01. From this expression it fol- 
27r; A02 = 



l^il 



^= (h Adl = 

J mi 



mi 
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(6) 



Let us now remark on what are the physical roots of 
this flux fractionalization. Further examining the solu- 
tion for the vortex (A0i = 27r, A02 = 0), we can write 
the vector potential as A = i^^^|A(r)| where r mea- 
sures the distance from the core and is a unit vector 
pointing along the core. The magnetic field is then given 
by |B| = i^(r|A|). The equation for the current J can 
then be rewritten as 
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dr 
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For such a vortex, the solution for the vector potential is 



l*il' 



1*1 




(7) 



Indeed the magnetic field vanishes exponentially 
away from the vortex core with the characteris- 
tic length scale given by the magnetic field pen- 
etration length A = 



*iP 
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,l*il' 
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Also, in contrast to the 



mi ^ \ y mi m2 

Abrikosov vortex ^ , besides having the fractionalization 
of magnetic flux our vortex also features the neutral vor- 
ticity. This, in particular, can be seen by substituting the 
solution (0) into (|l]). Then, at length scales larger than 
the magnetic field penetration length from the vortex 
core, we have the following expression for the energy den- 

'1/2 



sity: F = ^ 



l*il^ 
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1*21-^ 



1*1 



3 



1 

2m2 



-1/2 



4-5 



Thus the en- 



ergy per unit length of the vortex (A0i = 2n, A(p2 — 0) 
is divergent. This is due to the fact that such a topologi- 
cal configuration necessarily induces in a two-gap system 
the neutral superflow. Indeed the above expression for F 
is similar to the energy density of a vortex in a neutral 

system with effective stiffness 



1*1 



F„. 



1 l^lP |^2P 
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This shows transparently the physical origin of the pres- 
ence in the system of a massless neutral boson; a topo- 
logically nontrivial configuration (A(j)i = 27r, A(/)2 = 0), 
besides having a current in the condensate also nec- 
essarily induces current in the condensate 5*2 . Albeit in 
such a configuration there are no gradients of (j)2, how- 
ever the two condensates are not independent but are con- 
nected by the vector potential. The admixture of an oppo- 
sitely directed supercurrent of the condensate ^'2 leads 
to the situation when the two supercurrtents partially 
compensate the induced by each other magnetic field. 
This leads to the existence of the effective neutral super- 
flow in the system (which in particular should dominate 
long-range interaction between vortices). Moreover it is 
the fact that the two currents partially compensate in- 
duced by each other magnetic field, which leads to the 
fractionalization of magnetic flux. We stress that the vor- 
tex solutions (A(/)i — 2tt, A(p2 — 0) in this model, albeit 
being topologically stable, however can not form as an 
energetically preferred state in external field. In an ex- 
ternal field the system will form the composite vortices 
(A01 = 2tt, A02 = 2tt) described above. 

There is however an experimental set up which should 
allow one to directly observe (A(/>i = 2tt, A02 = 0) ex- 
citations. That is, in a two - gap superconducting film 
described by (0) there should occur thermal creation of 
pairs of vortices and antivortices with neutral vorticity 
and fractional magnetic flux. The fact that this type of 
vortex features both confined magnetic field and long- 
range interaction due to neutral superflow will allow: (i) 
detection of the Berezinskii-Kosterlitz-Thouless transi- 
tion in a system of these vortices by standard experi- 
mental techniques like flux-noise measurements or mea- 
surements in an applied dynamic magnetic field (like in 
a system of Abrikosov vortices (ii) such measure- 

ments straightforwardly allow one to extract data about 
type of interaction between a vortex and an antivortex 
in the system. In spite of being by definition sensitive 
only to magnetic flux carrying vortices these measure- 
ments should give exactly the picture of the BKT transi- 
tion like in a neutral U{1) system. That is, the transition 
should be a pure BKT transition but not a "would be" 
BKT transition like in a system of Abrikosov vortices. 



As it is well known there is no true BKT transition is a 
system of Abrikosov vortices: because their interaction 
is screened by the Meissner effect and a BKT transition 
is replaced by a crossover which is observable only when 
penetration length is sufficiently large p^ . Moreover 
the BKT transition in a system of these vortices should 
be observable even in a type-I system (both in the limits 
r] = and when rj is large, where one has sine- Gordon 
vortices interacting with a linear potential p^). 

In conclusion, we studied vortices in a two-gap super- 
conductor. We have shown that the system allows sev- 
eral types of vortices, such as vortices carrying an arbi- 
trary fraction of the magnetic flux quantum, which have 
no counterpart in ordinary one-gap superconductors or 
two-gap neutral superfluids. This study also has inter- 
disciplinary interest because similar models are highly 
relevant in particle physics (| jl^ . 

It is a great pleasure to thank L.D. Faddeev, D. 
Gorokhov, A. J. Niemi, G. Volovik, K. Zarembo, 
V. Cheianov, Y.M. Cho, O.Festin, S. Girvin and P. 
Svedlindh for interesting discussions and/or comments. 
This work has been supported by grant STINT IG2001- 
062, the Swedish Royal Academy of Science and Goran 
Gustafsson Stiftelse UU/KTH. The author also acknowl- 
edges NorFA mobility scholarship and Thordur Jonsson 
and Larus Thorlacius for hospitality during his stay at 
University of Iceland. 



[1] A.A. Abrikosov Sov. Phys. JETP 5 1174 (1957) 

[2] For a recent review and citations see G. E. Volovik 

Proc.Nat.Acad.Sci. 97 2431 (2000). 
[3] H. Suhl, B.T. Matthias, L.R. Walker, Phys.Rev.Lett.3, 552 

(1959) 

[4] A.J. Leggett, Prog, of Theor.Phys. 36, 901 (1966) 

[5] T.D. Lee, Phys.Rev. D8 1226, (1973) 

[6] see e.g. L. Shen et. al. Phys. Rev. Lett. 14 1025 (1965); G. 
Binnig et. al. Phys. Rev. Lett. 45, 1352 (1980) 

[7] see e.g. F. Bouquet et. al. Phys. Rev. Lett. 87 047001 
(2001) Amy Y. Liu et. al., Phys. Rev. Lett. 87 087005 
(2001) P. Szabo et. al. Phys. Rev. Lett. 87, 137005 (2001) 

[8] T. Yokoya et. al. Science 294 2518 (2001) 

[9] K. Moulopoulos and N. W. Ashcroft Phys. Rev. B 59 
12309 (1999); J. Oliva and N. W. Ashcroft Phys. Rev. B 
30 1326 (1984); N. W. Ashcroft J. Phys. A129 12 (2000). 
[10] E. Babaev, L.D. Faddeev, A.J. Niemi Phys. Rev. B65 
100512 (2002) 

[11] L. Faddeev, Quantisation of Solitons, preprint IAS Print- 
75-QS70, 1975; and in Relativity, Quanta and Cosmology 
vol. 1, M. Pantaleo, F. De Finis (eds.), Johnson Reprint, 
1979 

[12] Similar model has been found to be relevant for spin-triplet 
superconductors E. Babaev Phys. Rev. Lett. 88, 177002 
(2002) 



4 



[13] 



[14] 



[15] 

[16] 
[17] 



More details will b e present elsewhere. See also E. Babaev 
cond-mat/0201547. E. Babaev and D. Gorokhov in prepa- 
ration. 

The energy of the fractional flux vortices in a such a sys- 
tem will however be finite in a mesoscopic sample. Also, 
in a granular superconductor, lattice effects may rise mass 
of the neutral 0(2) boson which will make the energy per 
unit length of this defect finite. 

The author thanks G.E. Volovik for pointing out this cir- 
cumstance. 

P. Minnhagen, Rev. Mod. Phys. 59, 1001 (1987) 
We should also stress that, as it was shown by Pearl, 
in a one-gap superconducting film, vortices also in- 
teract due to existence of a "mushroom-like" mag- 



netic field outside film surface. This effect gives the 
vortices an 1/r interaction at a distance larger than 
[penetration length]^ /[film thickness] from the vortex 
core, while at a distance shorter than that the vortices 
interact logarithmically: J. Pearl, Appl. Phys. Lett. 5, 65 
(1964); and in Low Temperature Physics-LT 9, edited by 
J.G. Daunt et al. (Plenum, New York, 1965), p. 566. In 
contrast, the fractional flux vortices considered by us have 
logarithmic interaction at all length scales. 
[18] L. Faddeev, A.J. Niemi, Nature 387 (1997) 58; Phys. 
Rev. Lett. 82 (1999) 1624; Phys. Lett. B 525, 195 (2002) 
W. S. Bae , Y.M. Cho, S. W. Kimm Phys. Rev. D 65, 
025005 (2002) Y. M. Cho, H. W. Lee, D. G. Pak Phys. 
Lett. B 525, 347 (2002) 



5 



